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TS$DIVCD = D COMPLEX®16 / D COMPLEX®16 DIVISION R 16-SEP-1986 01: AX/VMS 4- ° 
a808 Breer ia8e T]:378o PATANES Bacte YOR OO ans1 O8 cl, 


1 «TITLE OTSSDIVCD = D COMPLEX*16 / D COMPLEX*16 DIVISION ROUTINE 
§ eIDENT /1-001/ 3 File: OTSDIVCD.MAR 
4; 
5 peeeeeeererererreeerereeeeeeeeeeneeeeeeeeeeeeeeeeeeTeKeeeeeeeeeeeeereeeeeress 
:* te 
09 : :* corre leny (c) 1978, 1980, 1982, 1984 BY ® 
0 8 :* DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. * 
000 19 = ALL RIGHTS RESERVED. * 
3 Pa 
00 11 ;* THIS SOFTWARE IS ea UNDER A LICENSE AND MAY BE USED AND sores * 
4 \§ 3* ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE * 
0 15 ;* INCLUSION OF THE ABOVE Spd NOTICE. THIS SOFTWARE OR A THE * 
00 14 ;* COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY * 
0000 15 ;* OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY * 
44 16 :* TRANSFERRED. * 
3 f 
0900 18 ;* THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE * 
0000 19 ;* AN SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT * 
4 $9 :* CORPORATION. * 
: © 
444 ¢ 3* DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS * 
000 3* SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. * 
0000 4 ;* 7 
0000 5 ;* « | 
0000 § PRRRRARERE AAA ARTA AAA TARA AAAAAAAA REAR AAAAARAATA AAA AREKe AKA ee ee ere tee Test | 
0000 ; 
0000 8 ; | 
44 9; FACILITY: MATH LIBRARY 
000 0 ;+ 
0000 1; ABSTRACT: 
0000 ¢ 3 : : 
0000 3 Perform D COMPLEX*16 division 
0000 4; 
0000 5 ;-- 
0000 «36: | 
0000 He 3 AUTHOR: 
0000 8 ; Steven B. Lionel, 12-July-1979 
0000 39; 
0000 40 ; MODIFIED BY: 
0000 41 ; 
9000 4¢ | 
0000 43: | 


QTRSDIVCD 


nist 


sort 


petal Baeomrcenetg oivis 16y 8 "OrSEb=18R2 94:39:82 FaNAnte. Saesorsoivcp.mar;1 
ig -SBTTL HISTORY ; Detailed Current Edit History 
is : Edit History 
3 


; 1-001 = Adapted from OTSS$DIVC version 1-003. SBL 12-July-1979 
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= D COMPLEX®16 / D COMPLEX®16 DIVISION R 16-SEP-1984 01:53:20 VAX/VMS Macro V0s- 
DECLARATIONS meat a TRE Si cA + 
3 LSBTTL DECLARATIONS 
5; 
¢ { INCLUDE FILES: 
09 $9 ; EXTERNAL SYMBOLS: 
& e 
2 86 5 
000 4 ; MACROS: 
00 66 ; 
00 67 ; 
000 68} PSECT DECLARATIONS: 
oooopesy «90 .PSECT _OTSSCODE PIC, USR, CON. REL. LCL, SHR, 
0000 3671 EXE. RD, NOWRT, 
9000 72 
44 "3 3 
000 7% : EQUATED SYMBOLS: 
00000006 9000 7a" 4 L part of dividend 
a = iy 
0000000¢ 0000 77 Ns = 12 | tmae pert of dividend 
00000014 0000 78 : = 20 eet bert of diviene 
0000001C 0000 79 d = 28 : imag part of divisor 
0000 80 
0000 «81: 
0000 82; OWN STORAGE: 
0000 «83 : 
0000 He : none 
0000 «8S 
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«SBTTL D COMPLEX*16 / D COMPLEX*16 giving D COMPLEX*16 result 


oT Esplvco = D COMPLEX*1 
1 D COMPLEX*16 


om 
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;+4 
; FUNCTIONAL DESCRIPTION: 
OTSSDIVCD_R3 = D COMPLEX*16 / D COMPLEX®*16 giving D COMPLEX*16 result 


The COMPLEX*16 result is computed as follows: 

1) Let (a, b) represent the COMPLEX*16 dividend. 

$} Let (c, d) represent the COMPLEX*16 divisor. 
) Let (r, i) represent the COMPLEX*16 quotient. 

Then: 


r = (ac + bd) / ¢ 
i = (be = ad) / (¢ 


CALLING SEQUENCE: 

Complex_quotient.wdc.w = OTSSDIVCD_R3(dividend.rdc.v, divisor.rdc.v) 
INPUT PARAMETERS: 

Dividend and divisor parameters are represented as 

FORTRAN D COMPLEX*16 numbers and are CALL BY VALUE. 

Passing 128 bit quantities by value is a violation 

of the VAX calling standard, but is excused because 

this is a code support routine not meant to be 

callable by users. 


IMPLICIT INPUTS: 
NONE 


DOOOOCOCOOOOOCOO 
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OUTPUT PARAMETERS: 
NONE 

IMPLICIT OUTPUTS: 
NONE 


FUNCTIONAL VALUE: 


The D COMPLEX*16 yplue returned is (a, b) / (c, d) 
in resisters RO-R3! This is a violation of the VAX 
calling standard, but is excused because this is 
: code support routine, not meant to be callable 

y users. 


SIDE EFFECTS: 


Modifies registers RO-R3! 

ither argument is a reserved operand. 
ogCing overflow 
Vv y 


ide zero 


FFP AAAI NIAIGEIPPIPIDINININPIDINID 2 9 SH QOOCOCOCOCOCOCOCOCOOOOOOOOOC@OOD ~o 


MO ODNAUE WIN 9 OD NAME WIN = OOO NAU EWN 9 OD NOAUE WR OOO NAOUS WO Oo~n 


0 
0 
0 
0 
0 
0 
0 
0 
0 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
0 
0 
0 
00 
00 
00 
00 
00 
00 


SOSCSSCSSSSSSSOSOSOSOSOSOSOOOOSOOSOSOSOSOSOSOSOOOSOOSOOSOOOSoOO 
a a ee a ed dd od 


SOOOOOoooooo 


nun 
II % 
wane 
~4D 
rem 

st 

< 

ns 

—_ 
~~ 


F 6 
OTSSHIVCD = D COMPLEX*16 / D COMPLEX*16 DIVISION R 19360-1386 91:33:80 AX/VMS Macro V04-00 Page ; 
1 D COMPLEX*16 / D COMPLEX*16 giving D COM 6-SEP-1984 11:27:54 CMTHRTL.SRCJOTSDIVCD.MAR;1 (5) | 
OF FO 144 LENTRY OTSSDIVCD_R3, “M<R4,R5,R6,R7,RB,/RO,R10,R11> | 
145 MTHSFLAG_ JACKET ; establish math error handler 
6D 00000000'GF 9€ MOVAB G*MTHSSJACKET_HND, (FP) | 
; set handler address to jacket 
9 3; handler 
3 14 
4 149 ; Perform scaling of all operands before division 
51 14 AC 8 7 EF 9 128 EXTZV a7, ae c(AP), R1 ; R1 = c(AP)<exp> eleetearr 
50 1C AC 8 7 «CF o 1 9 EXTZV #7, #8, d(AP), RO 3 RO = d(AP)<exp> alaseoat? 
1 9 81 ; 1 CMPW RO, R1 ; RO = MAX (c<exp>, d<exp>) 
14 1 1 ¢ BGTR 2$ 
28 1 0 1A 15 MOVW R1, RO 
0 ° E 1D 1546 2$: MNEGB 0. RO ; RO = scaling exponent 0,377,376,...,1 
50 50 9 4 ? 133 cor 4 RO, RO ; build a floating scale factor 
6 138 : scale all operands 
58 14AC 50 65 0026 138 MULD RO, c(AP), RB 3 RB-RO gets c 
an a a 086 180 MULDS RO. SCAPD. Ra” 7 es 
- a ° : a 
52 Oc AC 50 665 ; ; 193 MULD RO, DCAP), R2 33 
50 58 4 $3 Ben 188 MULD R4, RB, RO ; RO = ac 
56 SA C5265 164 MULD3 R2. R10, R6 : R? = bd 
0 6 60 004 165 ADDD R6, RO ; RO = ac+bd 
54 OSA 64 (004 166 MULD2 R16, R4 : R4 = ad 
3 38 64 0048 16 MULD2 R8, R ; Re = be 
4 62 48 168 SUBD R4, R 3; R2 = be - ad | 
5 58 64 ef 169 MULD R8, RB 3; RB = cc | 
5a SA O64 «(0051~Ss«170 MULD2 R10, R10 : R10 = dd | 
58 SA 60 0054 £4171 ADDD R10, RB ; RB = cc + dd 
50 58 66 0057 17¢ DIVD2 R8, RO ; RO = (ac+bd) / (cc+dd) 
52 58 66 2A V7 Divd » R2 : R2 = (bc-ad) / (cc+dd) 
5D 175 RET 3 (RO-R1, R2-R3) = (r, i) 
5E 178 
5E 17 END 
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MIHSSJACKET_HND feeeeeee 1 
OTSSDIVCO_R3 00000000 RG $1 
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! Psect synopsis ! 


beer esenrecconcaesad 


PSECT name Allocation PSECT No. Attributes 


- ABS , 94 44444 ( 0.) 89 ( 9} NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
-OTSSCODE OOOOOSE ( 94.) OTC 1.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 


tense sence w ese ween cen mere n + 


! Performance indicators ! 


beret eee etree eer eo oem mera > 


Phase Page faults CPU Time Elapsed Time 


Initialization 9 :00:00.0 0:00:04 

50 nee processing 121 09: -6 $: 
Synbol table sort 3 : : $9 3: a i 
cose f table output :00:00.01 0: a 
Psect synopsis output :00:00.0 0: 
Cross-reference output :00:00.0 0: 0009 
Assembler run totals 286 :00:01.9 ‘00; 11 


the working set Limit was 900 pages. 
3071 bytes (6 pages) of virtual memory were used to buffer the intermediate code. 
There were 10 pages of symbol table space oh ecenee to hold 6 non-local and 1 Local symbols. 
source Lines were read in Pass 1, producing 1 1 object records in Pass 2. 
page of virtual memory was used to “define 1 macro. 


eee weet eee mem mewn ene eee ewe ce e 


! Macro library statistics : 


Senensenensseennn cere woar eco s + 


Macro Library name Macros defined 
"$255$DUA28: [SYSLIBJSTARLET. MLB; 2 0 

0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESS1ON/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:OTSDIVCD/OBJ=OBJ$:OTSDIVCD MSRCS$:MTHJACKET/UPDATE=(ENHS:MTHJACKET) #MSRC 
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